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The theoretical status of Higgs boson production at the LHC within 
the minimal supersymmetric extension of the Standard Model is re- 
viewed. The focus will be on the evaluation of higher-order correc- 
tions to the production cross sections and their phenomenological 
implications. 
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The theoretical status of Higgs boson production at the LHC within the minimal 
supersymmetric extension of the Standard Model is reviewed. The focus will be 
on the evaluation of higher-order corrections to the production cross sections and 
their phenomenological implications. 

1 Introduction 

The search for Higgs bosons is one of the most important goals for present 
and future experiments. Once a Higgs will be found, its properties have to be 
investigated in order to distinguish between the Standard Model [SM] and its 
extensions. Supersymmetry represents one of the most attractive extensions, 
since it provides a natural solution to the hierarchy problem of the SM. The 
minimal supersymmetric extension of the SM [MSSM] contains two isospin 
doublets of Higgs fields, which are necessary to allo"«L to introduce up and 
down-quark masses without breaking supersymmetry El. Moreover, they are 
required to cancel anomalies associated with Higgsinos, the supersymmetric 
partners of the Higgs bosons. After electroweak symmetry breaking, one is left 
with a spectrum of 5 elementary Higgs particles: two neutral CP-even (/i, H), 
one neutral CP-odd {A) and two charged (H^) Higgs bosons. The MSSM 
Higgs sector can be described by two independent parameters at leading order 
(LO), which are in general chosen as tg/3 = V2/V1, i.e. the ratio of the two 
vacuum expectation values of the scalar Higgs fields, and the pseudoscalar 
mass Ma. Radiative corrections to the Higgs masses and couplings are largeH, 
since the leading part e grows as the fourth power of the top-quark mass rrit- 



where Gp denotes the Fermi constant and TOj the common squark mass. They 
shift the upper bound on the light scalar Higgs mass to < 130 GeV. The 
Higgs couplings to top (bottom) quarks and gauge bosons are modified by 
SUSY factors, which are collected in Table |^. The angle a denotes the mixing 
angle between the scalar Higgs particles h, H. An important property of the 
SUSY couplings is the enhancement (suppression) of the bottom (top) Yukawa 
coupling for increasing tg/3. The direct search for the MSSM Higgs particles 
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at LEP yields lower limits Mh^H,A ^ 62.5 GeV for the neutral Higgs massesH. 

Table 1: Higgs couplings in the MSSM to fermions and gauge bosons [V = W, Z] relative 
to SM couplings. 
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At the LHC, the production of neutral Higgs boson^ is dominated by 
gluon fusion gg — > [0 = /i, if, vl] □. Only for large values of tg/3 does Higgs 
bremsstrahlung off bottom quarks, gg, qq — > 4>bb, become dominant. This is 
shown in Fig. in which the Higgs boson production cross sections via the 
various mechanisms for the scalar and pseudoscalar Higgs particles at the LHC 
are presented. For their discovery, several decay jnodes must be exploited in 
different regions of the MSSM parameter space h — + 77; H,A ^ t^t^; 
H ^hh^ bb-fj; ZH ^ i+£-bb; H,A^ tt. 



2 gg-^ (j> 

2.1 Lowest order 

At LO the gluon fusion mechanism is mediated by heavy 
squark loops, see Fig. 0. The LO cross sections are given by 



<^Lo{pp 
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t Charged Higgs bosons will be produced via radiation from a top quark|i or in pairs via 
the Drell-Yan processa or gluon-gluon collisionsLl. They will not be considered here. 
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Figure 1: Neutral MSSM Higgs production cross sections at the LHC [-^/s = 14 TeV] 
for gluon fusion gg — > 0, vector-boson fusion qq — > qqVV — > qqh/qqH, vector-boson 
bremsstrahlung qq — > V* — > hV/HV and the associated production gg,qq — > cf>bb/4>tt in- 
cluding all known QCD corrections for tg/3 = 1.5,30. 



with the scahng variables — M^/s and 



4m^ -/Af ? and the function 



/(r) = 



LOS , — tTT 



r > 1 
T < 1 



(3) 



C^^ denotes the gluon luminosity. The top and bottom quarks as well as 
the third-generation squarks provide the dominant contributions to the cross 
sections. Because of the behavior of the SUSY couplings the top (bottom) con- 
tributions are suppressed (enhanced) for large tg/3. The squark contributions 
are sizeable only if the squark mass is below ~ 400 GeV, as can be seen in 
Fig. H, where the ratio of the light scalar Higgs production cross sections-ssath 
and without the squark loops is shown as a function of the squark masst2ltil. 
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Figure 2: Diagrams contributing to gg H at lowest order. 
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Figure 3: Ratio of the cross section cr{pp — » h + X) with and without squark loops as a 
function of the common squark mass rrig for two values of tgf3 = 1.5,30, and for Ma = 
100 GeV. The secondary axes present the corresponding light scalar Higgs mass Mf^. 

2.2 QCD corrections 




Figure 4: Typical diagrams contributing to the virtual and real QCD corrections to gg — > H. 

Quark loops. The NLO QCD corrections consist of two-loop virtual cor- 
rections and one-loop real corrections from the processes gg — + Hg,gq 
Hq, qq — > Hg. Typical diagrams are depicted in Fig. ^ The evaluation of 
the virtual and real corrections has been performed within dimensional regu- 
larization. The five-dimensional Feynman integrals of the two-loop diagrams 
have been reduced analytically to one-dimensional ones, which have been inte- 
grated numerically. The tensor reduction of the virtual three-point functions 
cannot be performed completely down to scalar integrals, but one is left with 
the calculation of irreducible tensor integrals. The heavy quark mass has been 



4 





K(pp ^ h/H + X) 
Vs= l4TeV 


H = M ^ M|^^ 
M, - 175 GeV . 
CTEQ4M 








m_ 


tgp-1.5 




m 




m 

® 


(H) 


m 



K(pp ^ A + X) 


H = M ^ 


Vs=14TeV 


M, = 175 GeV 




CTEQ4M 








!■-;■, 















Figure 5: K factors of the QCD-corrected gluon-fusion cross section cr(pp ij> + X) at the 
LHC with cm. energy yfs = 14 TeV. The dashed lines show the individual contributions of 
the four terms of the QCD corrections given in eq. M). The renormalization and factoriza- 
tion scales have been identified with the corresponding Higgs mass, fi = M = M^, and the 
CTEQ4M parton densities have been adopted. 

renormalized on-shell and the strong coupling as in the MS scheme. The re- 
maining coUincar singularities of the real corrections are absorbed in the NLO 
parton densities, defined in the MS scheme. In order to calculate the QCD 
corrections to the pseudoscalar Higgs boson production, a consistent scheme 
for the Jife coupling has to be used_We implemented the 't Hooft-Veltman 
scheme Ej, its modification by LarinEj and the scheme introduced by Kreimer 
li^, all schemes giving identical results; for a naive 75 the final result turned 
out to be ambiguous and thus inconsistent, as a result of the ABJ anomalyll3. 
The cross sections can be cast into the formu 
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Figure 6: The renormalization and factorization scale dependence of the Higgs production 
cross section at lowest and next-to-leading order for two different Higgs bosons H, A with 
masses Mjj = 500 GeV and Mj\ = 100 GeV and two values of tg/S = 1.5, 30. 
r-1 



dry -j^ ^ dqq{z,TQ} 



1 



with z = 'uJt — M^/s. Pgg and Pgq are the standard AltareUi-Parisi sphtting 
functions E^. The K factors K = (Tnlo/'^lo^ where the NLO (LO) cross 
sections have been evaluated with NLO (LO) as and parton densities, are 
presented in Fig. || as functions of the corresponding Higgs masses. The QCD 
correctious are large and positive, increasing the production cross sections by 
10-100%Q, so that they can no longer be neglected. In spite of the large size of 
the QCD corrections the scale dependeHce is reduced significantly, rendering 
these NLO results reliable within ~ 20%cl, as can be seen in Fig. ^, which shows 
the LO and NLO cross sections as a function of the common renormalization 
and factorization scales in units of the corresponding Higgs masses. Improved 
measurements of the parton densities in the HERA experiinents reduced the 
uncertainties from the parton densities to a level of ^ 15% 0. The size of the 
K factors depends strongly on tg/3B. 

In the limit of heavy quark masses the mass-dependent terms of Eq. ^ 
reduce to very simple expressionsQ: 



„h,H 



C^iTQ) ^ 6 



dgq{z,TQ) -Z^ - (1 - zf 



dgg{z,TQ) -y(l - zf 
32 

dqqiz,TQ)^ —{I-Zf (5) 



These limits can also be derived from low-energy theorems. 



Low-energy theorems. For scalar Higgs bosons these rely on the fact that for 
vanishing Higgs momentum the entire interaction of the Higgs particle with 
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fermions and gauge bosons can be generated by a shift of masses by the Higgs 
field. Thus matrix elements with an externajJjght Higgs boson are related to 
the matrix element without the Higgs bosonEj: 

lim MoiXH) = ^moJ-MoiX) (6) 

where X denotes any particle configuration, v — 1/ \/ \/2G f the vacuum ex- 
pectation value, and g'^ the corresponding SUSY coupling of Table |l} In order 
to extend this relation to higher orders in perturbation theory, it must be for- 
mulated in terms of bare quantitiesla. For on-shell Higgs particles this mathe- 
matical limit coincides with the massless Higgs limit. Thus we can derive the 
effective coupling of a light Higgs boson to gluons from the gluon self-energyu, 

£,„ = J|-G-G:„f{l + fiii + (7) 

This Lagrangian has to be interpreted as part of the basic Lagrangian describ- 
ing the effective theory, once the top quark is integrated out, and thus provides 
the proper description of the heavy quark limit. The effective coupling has to 
be inserted in the blobs of the effective diagrams shown in Fig. 0. The final 
result coincides with Eqs. of the scalar Higgs bosons. 




Figure 7: Typical effective diagrams contributing to the QCD corrections to gg H . 

The low-energy theorems for pseudoscalar Higgs particles are based on the 
ABJ anomaly in the divergence of the axial- vector current 113, 

d^jl = 2mQQz75Q + ^G^^^Gl, (8) 

ZTT 

with G"^ = ^e^iyct/jG""'' denoting the dual field strength tensor. The axial- 
vector current operator fulfils the low-energy condition EEI 

{gg\d^,3^A\A) for -> (9) 

Using the basic interaction Cint = ~S.Q''^QQ^lbQ -^1 ^ the effective Lagrangian 
for the Agg coupling can be derivedu: 

t^e}S=4^G'^'''Gl^^^ (10) 
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Due to the Adler-Bardeen theorem of the non-renormalization of the ABJ 
anomaly t3, this effective Lagrangian is vahd up to aU orders of perturbation 
theory. Its insertion into the effective diagrams analogous to Fig. ^ leads to 
the results of Eqs. of the pseudoscalar Higgs boson. 

Squark loops. If the full massive LO cross section is multiplied by the K factor 
obtained in the heavy quark limit, a reliable approximation is obtaineid- within 
10% for the top quark contribution to the production cross sections llil. Thus 
it will be sufficient to derive the QCD corrections to squark loops in the heavy 
squark limit. This can be done by extending the scalar low-energy theorems 
to squarks. We performed this calculation in the approximation of degenerate 
squark flavors and very heavy gluinos so that Ihe latter decouple. The effective 
Lagrangian mediated by a squark loop reads 

^eff = ^G^^'-G^.- |l + f ^ + 0{al)] (11) 
487r V [ on J 

The amplitudes can be obtained from effective diagrams in analogy to Fig. ^ 
They were added to the massive quark amplitudes and contracted with the 
massive LO form factors in the virtual corrections. In this way we derived the 
most reliable approximation. The K factors differ only by less than ^ 10% 
from the K factors for the quark loops alone iB. This is shown in Fig. |[ 
which presents the K factors with and without squark loops for a common 
squark mass of 200 GeV, for which the LO cross sections are significantly 
enhanced. Thus the full K factors can simply be approximated by the massive 
K factors to the quark loops alone, while the LO cross sections should include 
the massive squark contributions. This approximation is valid, because the 
QCD corrections for heavy particle loops are dominated by soft and coUinear 
gluon effects. 

2.S Soft gluon resummation 

The results of the last paragraph provide a strong motivation for the resum- 
mation of soft gluon effects in the Higgs boson production cross sections. This 
step has been carried out for a common squark mass TOq = 1 TeV, so that 
squark loops can be neglected 0. Moreover, we worked in the heavy quark 
limit, which provides an approximation to the cross sections within ~ 25% for 
tg/3 < 5. For larger values of tg/3 the bottom contribution becomes significant. 
Finally the gq and qq initial states will be neglected, since their contribution 
only amounts to ~ 10% Within this approximation the partonic cross sec- 
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Figure 8: K factors of the cross sections cr{pp h/H + X) with [solid lines] and without 
[dashed lines] squark loops as a function of the corresponding scalar Higgs mass for two 
values of tg(3 = 1.5,30. The common squark mass has been chosen as M~ = 200 GeV. 

tions factorize into three pieces0, 

^gg = <^qI^<PP4> (12) 
The factors are fixed by the effective Lagrangian at NNLO00: 

HA = 1 (13) 

lla^(Mt) 3866 - 201 A^F /a, (Mt)^ " 
Hh,H = 1 + 1- - 



TT 144 
153-197Vf a,{MH) - a.iMt) , 3, 



33 - 2Nf 



0{at) 



while the correction factors originate from effective diagrams. At LO they 
are normahzed as — (5(1 — z) with z — ]\4^/s. It has been shown for the 
DreU-Yan process that for z < 1 the cross section factorizes into soft gluon 
contributions, jet functions containing coUinear gluon effects and a hard matrix 
elementJ23. From this factorization the Sudakov evolution equation can be 
derivedc3: 

r2 



Mi 



dMl 



(14) 



where the convolution is defined ds, f ® g = dz' J{z')g{z/ z'). The evolution 
kernel can be determined from the perturbative result by matching the 
leading terms with the perturbative expansion of the solution to Eq. In 
this way we have determined the first term in the expansion of and thus 
resummed the soft gluon effects with NLO accuracy. The renormalization and 
mass factorization defined the strong coupling as and the parton densities in 
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the MS scheme. The resummed result has then been used as a generating 
functioiiaJ for the perturbative expansion up to NNLO. The NLO terms of 
read asE^I 



= 12I?i(z) - 24f i(z) - 6Vn{z)L^ + n'S{l - z) 
Pa^ = pi/if + 65(1 



(15) 
(16) 



where 



V,{z) = 



1- z 



f,(z) = log^(l-z), L,^logi^\ (17) 



The term £1(2;) extends the conventional resummation techniques, which only 
resum the soft gluon logarithms 'Di{z). The leading part of the £i{z) terms 
is of a pure coUinear nature and thus universal as well. The consistency of 
this extension has not been proved so far. However, it is important to note 
that these terms are large for LHC processes, and thus relevant. Moreover, 
the analogous analysis for the Drell-Yan process and deep inelastic scattering 
reproduced the leading terms in £i{z) at the NNLO level so that we got con- 
fidence that they are universally factoriziii&in the coUinear limit. The NNLO 
expansion of can be cast into the formcj 



= 3 (24^3(2) + (-2/5o - mL^)V2{z) + (-24C2 + 2/?oL^ + l2Ll)V^{z) 



(48C3 + 12(2^^ - -(3oLf,)Vo{z) - iSSsiz) 



+ (4/3o + 24 + 72L^)£2iz) + (48(2 - 4/3oi^ - 24L^ - 24Ll)£,{z) 
2909 

+ (I8C2' - 36C4 - ^/3o + C2/3oiM - 24(3^^ - 6C2i^)<5(l - z) 



Pa 



pf^^ + 3 {24I?i(z) - l2L^V^{z) - 48£i(z) 
+ (12C2 + 6 + /5oi^)^(l-z)} 



(18) 



with /3o = (33 — 2Np)/6. The analogous analy^ of the DrcU-Yan process 
results in a reliable approximation up to NNLOEj so that the NNLO results 
of Eq. m are expected to be a reliable approximation of the NNLO corrections 
to the Higgs production cross sections. The correction factors convoluted with 
NLO parton densities are presented in Fig. ^as a function of the corresponding 
Higgs mass. The NLO expansion of is denoted by 71 and the NNLO ex- 
pansion by 72. We observe a good agreement of 71 with the exact NLJD result, 
within ~ 5%, while the correction factor 72 turns out to be large This, 
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however, is caused by the use of NLO quantities in all perturbative orders of 
the correction factor. A consistent analysis requires LO and parton densi- 
ties for the LO correction factor and NNLO quantities for 72. This reduces the 
NLO correction factor to a level of ^ 1.5, which is significantly smaller than 
the NLO and 71 curves in Fig. Thus in order to obtain a NNLO prediction 
NNLO parton densities are needed, which, however, are not available so far. 
Therefore a consistent theoretical prediction of the NNLO Higgs production 
cross sections is not yet possible. 
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Figure 9: Exact and approximate two- and three-loop correction factor convoluted with NLO 
gluon densities in the heavy top quark limit for the MSSM Higgs bosons. The CTEQ4M 
parton densities have been adopted with as{Mz) = 0.116 at NLO. 



3 Conclusions 



We have reviewed the theoretical status of the production of neutral MSSM 
Higgs bosons at the LHC. The NLO QCD corrections are large and positive, 
while the theoretical uncertainty estimated from the residual scale dependence 
is reduced significantly. The squark loop contributions become sizeable for 
squark masses below 400 GeV. The K factors, however, deviate by less than 
~ 10% if squark loops are included. The K factors depend strongly on the 
MSSM parameter tg/3. Finally we have performed the soft gluon resummation 
for the Higgs production cross sections, which we expanded up to NNLO in 
order to get a reliable estimate of the NNLO corrections. These turn out to be 
potentially large. However, a consistent phenomenological analysis at NNLO 
requires NNLO parton densities, which are not yet available. 
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